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Abstract 

We construct a DGP-inspired braneworld model where induced gravity on the brane 
is modified in the spirit of f{R) gravity and stringy effects are taken into account by 
incorporation of the Gauss-Bonnet term in the bulk action. We explore cosmological 
dynamics of this model and we show that this scenario is a successful alternative for dark 
energy proposal. Interestingly, it realizes the phantom-like behavior without introduction 
of any phantom field on the brane and the effective equation of state parameter crosses 
the cosmological constant line naturally in the same way as observational data suggest. 
PACS: 04.50.-h, 04.50.Kd, 95.36.-hx 
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1 Introduction 



According to recent cosmological observations, our universe is undergoing an accelerating phase 
of expansion and transition to the accelerated phase has been occurred in the recent cosmologi- 
cal past [1] . In order to explain this remarkable behavior and despite the intuition that this can 
be achieved essentially only through a fundamental theory of nature, we can still propose some 
theoretical approaches to describe this astonishing feature. It is thought that the present accel- 
eration of the universe expansion might be due to a dark energy component in the universe [2]. 
The simplest way to describe the accelerated universe is to use the cosmological constant in the 
Einstein field equations. However, huge amount of fine-tuning required for its magnitude and 
other theoretical problems such as unknown origin and lake of dynamics make it unfavorable for 
this purpose. It is then tempting to find alternative models of dark energy which have dynami- 
cal nature unlike the cosmological constant. In this respect, it is a challenging job in theoretical 
physics to formulate a theory for cosmological evolution with capability to address the origin 
of dark energy too. It is known from cosmological observations that the dark energy content of 
the universe has a major contribution on the total energy budget of the universe. Nevertheless, 
the usual fields available in the standard model of particle physics are not sufficient to account 
for huge dark energy reservoir in the universe. In this respect, the concept of dark energy in 
the Einsteins gravity just with standard matter or fields cannot be implemented. Therefore, a 
modification to the Einsteins field equations either in the geometric or in the matter sector is 
essential to accommodate the present cosmological observations. In this viewpoint, Einstein's 
field equations can be reformulated as either G^,^ = 87tG{T^^^^ + T^f "'''^)) ( where T^^^ and 
T'^"-^^^ are energy- momentum tensor of ordinary matter and dark energy respectively), or the 
geometric part of these equations acquire a modification to incorporate the " Dark Geometry" 
as G^j, + G^^"'"*^^ = S7rGT^^\ Both of these approaches have attracted a lot of attention in the 
recent years. 

In this respect, braneworld scenarios belong to the second category: they modify essentially 
the geometric sector of the Einstein's field equations. Within a braneworld viewpoint, the 
Dvali-Gabadadze-Porrati (DGP) setup provides a new mechanism to explain the late-time 
acceleration of the universe based on a modification of gravitational theory in an induced gravity 
(IG) perspective [3]. This new mechanism to explain cosmological late-time acceleration is more 
interesting than introduction of an exotic dark energy and it is able to explain the origin of the 
dark energy proposal from a purely gravitational viewpoint. Despite this astonishing feature, 
the self-accelerating branch of the DGP scenario is unstable due to existence of ghosts [4]. 
However, it seems that extension of the DGP scenario to more generahzed frameworks such as 
the model presented in this paper have the potential to overcome the ghost problem due to 
their wider parameter spaces [5]. On the other hand, if we accept the pure DGP setup as our 
gravitational theory, we still have some problem in order to explain all histories of the universe 
and to match in a smooth way different phases of the universe evolution. Fortunately, it has 
been shown that the normal, non-self- accelerating branch of the DGP scenario has the potential 
to explain the phantom-like behavior without introducing any phantom fields on the brane [6] . 
Recently, it has been shown that incorporation of the modified gravity with f{R) term on the 
normal branch of the DGP setup has the potential to realize the late-time acceleration [7]. This 
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interesting achievement opens new windows on the phantom-hke prescription in DGP-inspired 
scenarios, one of which is presented in this paper. 

From another perspective, in a braneworld scenario the radiative corrections in the bulk 
lead to higher curvature terms in the action. At high energies, the Einstein-Hilbert action 
will acquire quantum corrections. The Gauss-Bonnet (GB) combination arises as the leading 
bulk correction in the case of the heterotic string theory [8]. This term leads to second- 
order gravitational field equations linear in the second derivatives in the bulk metric which is 
ghost free [9,10,11], the property of curvature invariant of the Gauss-Bonnet term. Inclusion 
of the Gauss-Bonnet term in the bulk action results in a variety of novel phenomena which 
certainly affects the cosmological dynamics of these generalized braneworld setup, although 
these corrections are smaller than the usual Einstein-Hilbert terms [12-15]. In the presence 
of the Gauss-Bonnet (GB) term in the bulk action and induced gravity (IG) on the brane 
(GBIG), there are different cosmological scenarios, even if there is no matter in the bulk [16]. 
In this paper, we generalize the previous studies to the case that induced gravity on the brane 
is modified in the spirit of f{R) gravity. We show that there are several interesting features 
which affect certainly the cosmological dynamics on the brane. Since the Gauss-Bonnet and 
induced gravity effects are related to the two extremes of the scenario (ultra-violet (UV) and 
infra-red (IR) limits), inclusion of stringy effects via Gauss-Bonnet term leads to a finite density 
big bang. This interesting feature has been explained in a fascinating manner by T-duality of 
string theory [17]. 

The motivation for incorporation of the modified induced gravity on the brane lies in the 
fact that modified gravity emerges as serious candidate which is able to addressing the definitive 
answers to several fundamental questions about dark energy. For example, the origin of dark 
energy may be explained by terms that could be relevant at late times. Also these terms can 
be considered as the source of early time infiation [18]. Therefore, modified gravity is a natural 
scenario to have a unified theory to explain both the infiationary paradigm and dark energy 
problem. The standard Einsteins gravity may be modified at low curvature by including the 
terms that are important precisely at low curvature. The simplest possibility is to consider a 
term in the Einsteins- Hilbert action [19]. It has been suggested that such a theory may be 
suitable to derive cosmological models with late accelerating phase. Although a theory with 
-^-term in the Einsteins gravity accounts satisfactorily the present acceleration of the universe, 
it is realized that inclusion of such terms in the Einstein-Hilbert action leads to instabilities 
[20]. A modified theory of gravity which contains both positive and negative powers of the 
curvature scalar R namely, f{R) = R + aR"^ + (3^ where a and P represent coupling constants 
with arbitrary constants m and n are considered for exploring cosmological models [18]. It 
is known that the term R^ dominates and it permits power law infiation if 1 < m < 2, in 
the large curvature limit. Recently a phenomenologically more reliable form of f{R) has been 
proposed in Ref. [21] which we will consider as a suitable ansatz in our study. 

With these preliminaries, the motivation of this work is to construct a braneworld scenario 
with induced gravity and Gauss-Bonnet effect where the induced gravity itself is modified in 
the spirit oi f{R) gravity. Then we examine cosmological dynamics in this setup and we show 
that its account for phantom-like behavior without introducing any phantom field that violates 
the null energy condition neither on the brane nor in the bulk. This model naturally realizes a 
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smooth phantom divide hne crossing by its equation of state parameter and this crossing occurs 
in the way that is supported by observations, that is, from above —1 to its below. We show 
also that the phantom-like behavior occurs in this setup without violation of the null energy 
condition at least in some subspaces of the model parameter space. 



2 The Setup 

The action of our modified GBIG model contains the Gauss-Bonnet term in the bulk and 
modified induced gravity term on the brane as follows 



J brane 



M2 



(1) 



where (3{> 0) is the GB coupling constant, rc(> 0) is the IG cross-over scale, k.^ is the five 
dimensional gravitational constant, A is the brane tension and A5 is the bulk cosmological 
constant. Also K is the trace of the mean extrinsic curvature of the brane defined as follows 



K, 



HI/ 



y=-e 
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(2) 



and its presence guarantees the correct matching conditions across the brane. We denote the 
matter field Lagrangian by with the following energy-momentum tensor 



The field equations resulting from the action (1) are given as follows 



(3) 



(4) 



The corrections to the Einstein Field Equations originating in the GB term are represented by 
the Lovelock tensor [22] 



A r>A -qK a -qKL -qA I 9 tdAKLM p ^ rA T 
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The energy-momentum tensor localized on the brane is 

M2 
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r't = -Mlf{R)G% - Rf{R) - f{R) + 2-^ 5^ + T^. 
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where a prime marks differentiation with respect to the argument and is the covariant 
derivative with respect to g^^^. The corresponding junction conditions relating quantities on 
the brane are as follows [23] 
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To formulate cosmological dynamics on the brane, we assume the following line element 

ds^ = t)dt^ + a^(y, t)-fijdx'dx^ + b'^{y, t)dy^ (8) 

where 'jij is a maximally symmetric 3-dimensional metric defined as = 6ij + k ^^^2 where 
k — —1, 0, +1 parameterizes the spatial curvature and = XiX\ y is the coordinate of extra 
dimension and the brane is located at y = 0. The junction conditions on the brane now gives 
the following expressions 



lim 



and 



dyU 



y=-e 



. _ 2nMl 



na 



k 



lim 



dya 



s^^' Mi 



fiR) 



+ 



n 



J y=0 



/'(i?)(2p(*''*)+3p(*°*)) 



y=0 



f'{R) 



3M| 



y=0 

(9) 



(10) 



If we choose a Gaussian normal coordinate system so that b'^{y,t) = 1, these equations with 
non- vanishing components of the Einsteins tensor in the bulk yield the following generalization 
of the Priedmann equation for cosmological dynamics on the brane 
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Rf{R) - f{R) - 6HRf"{R)^ 
where is the density of ordinary matter on the brane. $ is determined by [16] 

$ + 2/3$2 



6 a4 



(12) 



where C is the bulk black hole mass. In our forthcoming arguments we set C = and by 
A5 = 0, we find $ = and ^ = For simplicity, we choose $ = 0. Now, if we define the 

2 2 

cosmological parameters as Urn = -ru^, = T^jp^ = ^chHq and flcurv = ^ ■T'^ , then 



the Friedmann equation on the brane can be expressed as follows 



1 

IW) 
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where E{z) = and we have set A = 0. Evaluating the Priedmann equation at 2; = gives a 
constraint equation on the cosmological parameters of the model as follows 



1 -|- 2,y^j.^(\ -\- ~l~ ^curv • (■^^) 

We note that a contribution to this equation from f'[R)\^^Q has been normalized to unity. This 
constraint equation implies that the subspace with Qrn + ^curv < 1 is unphysical in the model 
parameter space. It is important to note that this model can be constraint in confrontation 
with observational data. As we will show this model mimics the ACDM model by setting 
Qrn ^ 0.26, Qcurv ^ 0.75, ^ 0.0001 and Qa ^ 0.01 . 

In this framework, the total energy density can be defined as follows 

Ptot = Pm + P curvi 

(15) 

where by definition 

Pcurv = ^ (RfiR) -f{R)- 6HRf"{R)) . (16) 
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The continuity equation for Pcurv is [24] 



?,HlVtMRf"{R) _3 



Pcurv ~l~ 3-ff(Pcurv ~^ Pcurv) r /■;/7-,\no ^ ) 

[f {R)\ 

where a dot denotes ^ and a prime marks Also Qm is the present day matter density 
parameter. So we can deduce 

Pcurv = —^\RRf\R) - 6HRf"{R) - 6HRf"{R) - mR^f"{R) 

+3RHf'{R) - 3Hf{R) - ISH'Rf'iR)^ + (18) 

Using equations (16), (17) and (18) we find the curvature fluid equation of state parameter as 
follows 

1 RRf"{R) - 6HRf"iR) - 6HRf"{R) - 6HR^f"'{R) 
3H Rf'{R) - f{R) - 6HRf"{R) 

2 H',nMRnR)a-' .^g. 

M^Hif'iRW RfiR) - f{R) - QHRf"{R) ■ ^ ' 
On the other hand we know that 

p^ + 3//(p^+pJ = 0, (20) 

where 

Pm = POm{ — f = POm{l + zf . (21) 

d 
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So the total equation of state parameter is given as follows 

1 RRf'iR) - mRf\R) - 6HRf"{R) - QHR^f"\R) - Qp^{^){a - Ha) 
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'-^ total — — 1 — 
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which includes the matter energy density and pressure too. Albeit, if we consider a braneworld 
without ordinary matter content, then the right hand side of equation (17) vanishes and we 
find 

^ 1 RRfjR) -GHRfiR) -6HRf"{R) -6HR'r{R) 

Rf'{R)- f{R)-6HRf"{R) 

According to recent observational data including the type la supernovae Gold dataset, there 
exists the possibility that the effective equation of state parameter evolves from larger than 
— 1 (non- phantom phase) to less than -1 (phantom phase, in which super-acceleration is real- 
ized) . This means that Weff crosses —1 line (the phantom divide line) smoothly. A number 
of attempts to realize the crossing of the phantom divide have been made in the framework of 
general relativity. For instance, we could mention scalar-tensor theories with the nonminimal 
gravitational coupling between a scalar field and the scalar curvature or that between a scalar 
field and the Gauss-Bonnet term, one scalar field model with nonlinear kinetic terms or a non- 
linear higher-derivative one, phantom coupled to dark matter with an appropriate coupling, 
the thermodynamical inhomogeneous dark energy model, multiple kinetic k-essence, multi-field 
models (two scalar fields model, quintom consisting of phantom and canonical scalar fields), 
and the description of those models through the Parameterized Post-Friedmann approach, or 
a classical Dirac field or string-inspired models, non-local gravity, a model in loop quantum 
cosmology and a general consideration of the crossing of the phantom divide ( see for instance 
[25] and references therein) and very recently crossing with Lorentz invariance violating fields 
[26] . In addition, there are interesting models of modified gravity realizing the crossing of 
the phantom divide too [27]. In the present paper, we study possibility of realization of the 
phantom-like behavior and phantom divide line crossing in our modified GBIG model without 
introducing any phantom matter on the brane or in the bulk. For this purpose, we adopt some 
reliable ansatz for f{R) and a{t) as follows. 



2.1 Phantom-Like behavior and crossing of the phantom divide 

By phantom-like behavior one means an effective energy density which is positive and grows 
with time and its equation of state parameter stays always less than —1. In this subsection we 
study possible realization of the phantom-like behavior in this UV/IR- complete theory. To do 
this end, we define effective energy density and effective equation of state parameter as basic 
ingredients of our analysis. To find the effective energy density we can use 

^' = y(Pm + Pe//) (24) 
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and then we rewrite our Priedmann equation (11) in this way to find 



Peff 



1 + ^OlH'^)H - Y(pm + Pcurve) 



rKlf'{R) 
Using the continuity equation 

Peff + 3//(l + UJeff)Peff = , 



- Pv 



(25) 



(26) 



since 
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(27) 



we find 



3H ' 



(1 + SaH^)H - 
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Hpv 



(28) 



Now to proceed further, we should specify the form of f{R). Extended theories of gravity based 
on four dimensional f{R) scenarios should follow closely the expansion of a ACDM universe 
[7,28,29] and also could have distinctive signatures on the large scale structure of the universe 
[30,31]. During the last few years, several methods have been proposed to reconstruct the shape 
of f{R) from observations [32,24]. This has been done, for example, by using the dependence of 
the Hubble parameter with redshift which can be retrieved from astrophysical observations [24] . 
Among these attempts, the models presented in Ref. [28] are used in our forthcoming arguments 
to study phantom-like behavior and crossing of the phantom divide line in our modified GBIG 
model. We proceed as follows: we start with the following observationally suitable ansatz for 

/(fl) = fl-fl. ,°'^/^£„ , (29) 



where 



^ ^'^l + (i?/i?e)"' 

R = -6(H + 2H^) 



(30) 

for an spatially flat FRW type universe. Both a and Rc are free positive parameters. We also 
adopt the following observationally reliable ansatz [33] 



V 



1 

1-1/ 



(31) 



and we translate all of our cosmological dynamics equations in terms of red-shift using the 
following relation between a and z 

l + z=^. (32) 

a 



Figure 1 shows the behavior of the effective energy density, Pe//, versus redshift. The effective 
energy density grows with time and pe// > always. These are necessary conditions for 
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phantom-like behavior but not sufficient: the effective equation of state parameter should stay 
below —1. Figure 2 shows the behavior of 1 + cUg// versus redshift. The universe enters the 
phantom phase in the near past and currently it is in the phantom phase. The transition to the 
phantom phase has occurred at ^ ~ 0.25. We note that based on some observations, crossing 
of the phantom divide line by the equation of state parameter occurs at z ~ 0.25. For instance 
the Gold sample mildly favors a crossing of the phantom divide line at 2; ~ 0.25 while no such 
trend appears for the SNLS data ( see for instance the paper by Nesseris and Perivolaropoulos 
inRef. [25]). 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 1: Variation of the effective energy density versus the redshift. 
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Up to now we have shown that this modified GBIG scenario reahzes the phantom-hke 
behavior without introducing any phantom matter on the brane or in the bulk. This phantom- 
hke behavior is a mimicry of ACDM scenario [6,7]. Now, the deceleration parameter defined 
as 

<l = -[§-2+^] (33) 

takes the following form in our model 



^ TRf'\R) - J.(p 



1 + 8aif2 + 2rHf'{R) 



(34) 



Figure 3 shows behavior of q versus z. The universe has entered an accelerating phase in the 
past at z ~ 0.84. The variation of H with redshift gives another part of important information 
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Figure 3: Variation of the deceleration parameter versus redshift. 
about the cosmology of this model. It is given by 

H _ 1 rRH'f"{R)-^{pm+Pcurv) , . 

Hi Hi l + 8aH^ + 2rHf'{R) ' ^ ^ 

Figure 4 shows the variation of versus redshift. Since H < always, the model universe 
described here will not acquire super-acceleration and big-rip singularity in the future. Since 
we have not introduced any phantom field on the brane or in the bulk, we expect that the 
null energy condition to be respected at least in some subspaces of the model parameter space. 
Figure 5 shows the status of the null energy condition in this setup. As this figure shows, 
phantom-like behavior occurs in this setup without violation of the null energy condition. We 
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Figure 4: Variation of versus the redshift. 
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z 

Figure 5: Variation of pe// + Pefj versus the redshift. Null energy condition is respected in the 
phantom phase of the model. 
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note that this model mimics a ACDM model in several respects and in this sense it is important 
to confront its parameter space with recent observation to see its viability. 

Finally we focus on the appropriate ranges of n that crossing of the phantom divide line 
can be realized. Our analysis shows that Ucurv as defined by (19) crosses the phantom divide 
line smoothly if 2.42 < n < 4.68. Also, ujtotai crosses the phantom divide line if 2.44 < n < 
4.74. And finally, uj as defined in equation (23), crosses the phantom line if 2.43 < n < 
5.15. In our model, this crossing occurs at z ~ 0.27 which lies well in the vicinity of the 
observationally supported value of 2; ^ 0.25. Figure 6 shows the behavior of these equation 
of state parameters versus z for n = 4. As this figure shows, crossing of the phantom divide 
evolves from quintessence to phantom phase in the same way as observations suggest. 

0.6 -] 1 

0.4 
0.2 
0.0 
-0.2 

(O 

-0.4 
-0.6 
-0.8 

-1.0 
-1.2 



Figure 6: Smooth crossing of the phantom divide line. 



3 Summary 

In this paper we have constructed a modified GBIG scenario: it contains Gauss-Bonnet term 
as the UV sector of the theory and Induced Gravity effect as IR side of the model. The induced 
gravity on the brane is modified in the spirit of f{R) gravity. The cosmological dynamics in this 
braneworld setup mimics a ACDM model in several respects: it realizes phantom-like behavior 
without introducing any phantom matter on the brane. The effective energy density increases 
with cosmic time and effective equation of state parameter crosses the phantom divide line 
smoothly in the same way as observations suggest. This crossing behavior occurs at z ~ 0.27 
which lies well in the vicinity of the observationally supported value of z ~ 0.25. The phantom- 
like behavior realized in this model happens without violation of the null energy condition in 
the phantom phase. Finally, to have smooth crossing of the phantom divide line by ujejf, the 
value of parameter n in the Hu-Sawicki model should lie in the range 2.42 < n < 4.68. 
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